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Abstract. We study the prospect of explaining the out- 
bursts of Soft X-ray Transients (SXTs) by the thermal- 
viscous instability in a thin disc. This instability is linked 
to hydrogen ionization and is significantly changed when 
irradiation of the disc by X-rays from the inner regions is 
included. We present the first numerically reliable, physi- 
cally consistent calculations of the outburst cycles which 
include the effects of accretion disc irradiation. The decay 
from outburst is governed by irradiation, as pointed out by 
King & Ritter (1998), leading to slow exponential decays. 
At the end of the outburst, the disc is severely depleted, 
which lengthens the time needed to rebuild mass to the 
critical density for an outburst. Despite this, the long re- 
currence times and quiescent X-ray luminosities of SXTs 
still require the inner disc to be replaced by another type 
of flow in quiescence, presumably a hot advection dom- 
inated accretion flow (ADAF). We include the effects of 
such truncation and show that the resulting lightcurves 
are conclusively similar to those of SXTs like A0620-00. 
We conclude that the two-a disc instability model pro- 
vides an adequate description of the outbursts of SXTs 
when both truncation and irradiation are included. The 
values for the viscosities in outburst and in quiescence are 
comparable to those used in CVs. We discuss the model 
in the context of present observations. 

Key words: Accretion, accretion discs. Instabilities, bi- 
naries: close, X-rays: general 

1. Introduction 

Low mass X-ray binaries (LMXBs) are systems in which 
a low-mass donor star loses matter to a neutron star or 
a black hole via Roche lobe overflow. The infalling gas 
forms an accretion disc where it diffuses towards the cen- 
tral body by losing its angular momentum through viscous 
interactions. Soft X-ray Transients (SXTs), a sub-class of 
LMXBs, show large amplitude X-ray and optical outbursts 
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during which their luminosity increases by several orders 
of magnitude within a few days. The peak is sometimes 
followed by an exponential decay on a timescale of a month 
before the system returns to its very low luminosity quies- 
cent state. This 'fast-rise exponential-decay' (FRED) be- 
haviour is considered to be typical of SXT lightcurves and 
is at the focus of our modeling effort. One should keep in 
mind, however, that many if not most of the SXT out- 
bursts do not have this shape. The recurrence time be- 
tween outbursts is from several months to tens of years 
(se e the reviews of Tanaka & Shibazaki 1996 and Chen et 
al. 1997 ). All of the known black hole LMXBs are transient 
systems. 

LMXBs and cataclysmic variables (CVs) show many 
similar characteristics and the analogy has led to the in- 
terpretation of the outbursts of SXTs as being due to a 
thermal-viscous instability in a geometrically thin, opti- 
cally thick Keplerian accretion disc, just as in dwarf no- 
vae (e.g. Cannizzo et al. 1982). In this disc instability 
model (DIM), a steady-state flow is stable if hydrogen is 
everywhere ionized. This would be the case for persistent 
LMXBs and their CV cousins, the novae-like systems. But 
if the temperature or (equivalently) the mass accretion 
rate from the donor star Mtr is low enough for hydrogen to 
recombine, the disc becomes thermally and viscously un- 
stable, oscillating between a hot, ionized state (outburst) 
and a cold, neutral state (quiescence). This is due to the 
strong dependence of the opacity on temperature when 
hydrogen is partially ionized (for reviews of the DIM, see 
Osaki |1996| ; Cannizzo |1998a| ; Lasota |200l"l ). 

Unfortunately, in this standard form, the DIM fails to 
provide an adequate model of LMXBs. For instance, the 
DIM predicts almost all LMXBs to be transient systems, 
in clear contradiction with observations. Van Paradijs 
( 1996| ) noted that the X-ray flux coming from (or close 
to) the compact object could heat the outer accretion disc 
above the hydrogen ionization temperature and therefore 
stabilize the flow. There is strong observational evidence 
of such irradiation e.g. the optical emission of LMXBs 
which comes mainly from X-rays reprocessed in the outer 
disc (van Paradijs & McClintock 1994). Yet, the simple 
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model consisting of a point source irradiating a thin disc 
does not lead to any significant heating of the outer disc 



(Dubus et al. 1999; see also Meyer & Meyer-Hofmeister 
1982| ). The geometry is certainly more complex and may 



involve indirect rather than direct irradiation of the disc. 
A corona reflecting only a small fraction 0.5%) of the 
X-ray luminosity from the center towards the outer disc 
would result in enough heating to quench the instability. 
Similarly, a disc warped out of the orbital plane could in- 
tercept enough of the central X-ray flux. 

Despite this gross uncertainty on the geometry of ir- 
radiation heating, there is no doubt that it plays a major 
role in the physics of LMXBs. The obvious next step in 
finding a model for SXTs, which we address here, is to 
examine the consequences of including irradiation heat- 
ing. This is particularly interesting since the standard 
DIM completely fails to explain the outbursts of SXTs. 
Numerical models show outbursts lasting at most only a 
few tens of days and with short recurrence times (Mi- 
neshige & Wheeler 1989; Menou et al. 2000). Moreover, 
the standard DIM predicts a rapid decay from the out- 
burst peak instead of a slower exponential decline. This 
was interpreted as a failure of the standard 'two-a' as- 
sumption (see §2.4 below) used to describe the viscous in- 
teractions in the DIM (Cannizzo et al. 1995). But King & 
Ritter (1998) pointed out that including irradiation heat- 
ing might naturally lead to the observed exponential de- 
cays without resorting to more complex formulations of 
the viscosity. Previous studies which argued that irradi- 
ation should have a negligible influence on the outburst 



properties (e.g. Cannizzo et al. 199?;, Cannizzo 1998b) as- 
sumed point source irradiation for which self-shadowing is 
crucial (Meyer & Meyer-Hofmeister 1982 ; Mineshige et al. 
1990| ). We have shown such arguments to be incorrect for 



persistent LMXBs and , in a ll likelihood, this also applies 
to SXTs (Dubus et al. |1999[). 



Furthermore, it is quite clear that the thin disc cannot 
extend down to the surface of the compact object in quies- 
cence: the X-ray detection of quiescent SXTs implies mass 
transfer rates onto the compact object which are four or- 
ders of magnitude higher than those predicted by the stan- 
dard DIM (Lasota 1996| ). This difficulty is resolved by a 
two-component accretion flow in which the inner thin disc 
gradually evaporates into an advection-dominated accre- 



tion flow (ADAF) during the decline (Lasota et al. 1996 
; for reviews, see Lasota 1999; Narayan, Mahadevan & 
Quataert 1999 ; note that a different model was suggested 
for neutron-star SXTs by Brown et al. 1998). This idea 



The outburst lightcurves predicted by the truncated 
disc instability models (TDIM) have been recently inves- 
tigated by Menou et al. ( 2000D ; (se e also Cannizzo, 1998b; 
Meyer & Meyer-Hofmeister, 1999] ). Menou et al. showed 
that, if the quiescent viscosity is smaller (by factors ~ 4) 
than the value usually assumed in dwarf novae, the TDIM 
results are in reasonable agreement with some of the ob- 
served properties of SXTs. Standard values of viscosity 
in quiescence lead t o ou tbursts which are too short, with 
reflares (see below §3T), and to short recurrence times. 
The major drawback of this model is, however, that the 
observed disc irradiation is not taken into account. 

In this work we show that a modified DIM including 
the effects of evaporation and irradiation can reproduce 
reasonably well the outbursts of SXTs with standard val- 
ues of the viscosity. The assumptions and the numerical 
technique used in the model are described in §2. The pre- 
dictions of the classical DIM are briefly recalled in §3. We 
successively explore the effects of irradiation heating with- 
out evaporation in §4 and with evaporation in §5. In §6 
we complete our investigation by varying different param- 
eters of the full model before discussing the results of this 
study in §7. 



2. Numerical model 



The numerical scheme used is that of Hameury et al. 
(1998), modified to include irradiation heating as in Dubus 
et al. (1999) and Hameury et al. (1999). Evaporation of 
the inner disc is treated in the same way as Menou et al. 

( poool ). 



Basically, the model takes advantage of the fact that in 
a thin disc one can decouple the vertical structure from the 
radial evolution equations. The vertical structure equa- 
tions are solved flrst, providing a grid from which the heat- 
ing and cooling terms can be obtained during the actual 
model computation. The radial disc equations are then 
solved on an adaptive mesh where the number of grid 
points is highest in the regions of steep temperature or 
density gradients. This results in great accuracy and elim- 
inates pure numerical uncertainties which plagued earlier 
attempts (e.g. Mineshige & Wheeler 198£ ; Cannizzo et al. 
1995 ). We discuss here some aspects related to the time- 
dependent irradiated disc model that were not discussed 
in the above references. 



2.1. Vertical structure 



has provided a good framework to understand the spectra 

of SXTs (e.g. Esin et al. 1997), their quiescent luminos- The series of disc vertical structures, including irradiation. 



ity (e.g. Narayan et al. 1996) and the X-ray/optical delay is calculated by solving eq. (22-27) of Dubus et al. (1999) 



during the rise of the outburst (Hameury et al. 1997). 
The evaporated disc does not participate in the thermal- 
viscous instability and, from a dynamical viewpoint, the 
disc can essentially be considered as truncated inside. 



assuming thermal balance and given values for the radius 
R, the surface column density E, the midplane temper- 
ature Tc and the irradiation temperature Ti„. In these 
equations, X-ray heating is limited to a thin layer above 
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the disc and appears in the boundary condition at the 
surface: 

^surf — "^cff + ^irr (1) 

Viscous heating follows the usual a prescription i.e. is pro- 
portional to_the local pressure in the layer P (Shakura & 
Sunyaev, 1973 ): 

^^i^foOKP (2) 

where fix = [GM\I 'B?f^'^ is the Keplerian frequency. 

The solutions give for each set of parameters unique 
values for the effective temperature Toff and a which are 
then stored. We compute 70 x 120 x 190 x 9 vertical 
structures in fix, S, Tc, Ti„ (Mi and R enter only via the 
Keplerian frequency) . We typically consider parameters in 
the range 10^-10^^ cm for R, 10^^-10^ g cm"^ ^j. iq3_ 
10^ K for Tc, and 0-Tc for Tirr. In most cases, this samples 
very well the range of values needed by the disc evolution 
code. 

For given Mi, R and Tirr, the solutions obtained above 
show an S shape when plotted in the S-Tc plane. Exam- 
ples of such irradiated S'-curves can be found in Tuchman 
et al. ( |1990D and Dubus et al. ( [19991 ). The low, cold branch 
corresponds to a layer of neutral hydrogen while in the 
upper, hot branch hydrogen is ionized. In between, the 
gradual ionization of hydrogen is associated with a steep 
increase of the opacity. This results in the middle branch of 
the S'-curve being thermally and viscously unstable. A ring 
of matter on the unstable branch will evolve either to the 
hot upper branch or the cold lower branch which are sta- 
ble. Neighbouring rings are affected somewhat like domi- 
nos and the instability propagates, periodically switching 
the disc between the hot state and cold state. This is the 
essence of the disc instability model (DIM) whose appli- 
cation to SXTs we study here. 

Following are useful numerical fits to the turning 
points of the S'-curve: Smax, Smin and the associated Tc 
when irradiation is taken into account. Emax is the max- 
imum density on the cold, neutral, branch while Smin is 
the minimum density of the hot, ionized, branch. We as- 
sume that a takes a different value on the hot (ah) and 
on the cold branch (etc ; see below §2.4). 

Smax = (10.8 - 10.30 

X a.-" «4Afr°-''+"-'^<""° ''^ gcm-2(3) 

-0.1 D-O-055 



Tc(S„,ax) = 10700 

^min ( 
-^c(^min) 



10 



K 



3 - 7.10 ah° ''^r° ''^' ''~° ''^ 



"-10 



g cm 



(4) 
' (5) 



(20900 - 11300O 



|1984 ). Irradiation thus has a much stronger influence on 

^^^(Sinin) than on Tc(Smax)- 

The S-curve disappears when Tc(Emax) > ?c(Sniin)- 
Irradiation heating is then enough to fully ionize the layer. 
This happens for Ti„ « 10^ K with a weak dependence on 
Ml, a and R (Tuchman et al. |1990[ Dubus et al. [19991 ). 
Note that this does not mean t he ve rtical structure be- 
comes isothermal (Dubus et al. 19991 ). A disc where Tirr 
is everywhere greater than 10* K will be thermally and 
viscously stable. 

2.2. Disc evolution equations 

The disc evolution is governed by the usual vertically inte- 
grated radial equations for mass, angular momentum and 

[l998| 



see 



energy conser vation : (Eq. 1-6 of Hameury et al 

also Cannizzo [1993[ ) to which we add an equation for Tj 
(Eq. [s] below, i j2.3| ). In particular, the energy conservation 
equation is: 

dTc ^ 2(Q+ -Q- + J) jRTe 1 djrv,) dTc 
dt CpY. fiCp r dr 9r ' ^ ' 

where (5+ and Q~ are the surface heating and cooling 
rates respectively. These are determined using the grid 
of vertical structures calculated above. The details of the 



procedure may be found in §3 of Hameury et al. 1998 



2.3. Irradiation temperature 

We have shown in Dubus et al. ( [1999[ ) that the simple 
assumption of a central irradiating point source in the 
plane of a thin disc led to results which were incompatible 
with observations. Following previous work (e.g. Shakura 
& Sunyaev 1973) we write the irradiation flux as: 



ctT 



C- 



with Lx = e min (Mjn, AfEdd) c 



(8) 



^ „-0.22 ,,r-0.01 pO.05-0.125 TA 

X a, M. R,n K 



""1 ^"-lo (6) 
where ^ = (Tirr/lO* K)^, Mi is in solar masses and 
i?io = i?/10^°cm. These fits are very close to those of 
Hameury et al. ( 1998[) in the case of no irradiation. Eq.(29) 
of Dubus et al. (1999) can be considered as a special case 
of Eqjl above for steady-state and a given dependence of 
Tiir on R. At the Srnin turning point, the maximum of the Mus) (Esin et al. 2000a; 2000b). Radiation-driven warp- 
opacity is reached close to the photosphere while it is at 
the bottom of the layer at Smax (Cannizzo & Wheeler, 



47ri?2 

Min is the mass accretion rate at the inner disc radius and 
-^^^Edd is the Eddington mass accretion rate for an accre- 
tion efficiency e of 0.1: MEdd = 1-4 x lO^^Mi gs'^ The 
maximum luminosity that can irradiate the disc is thus 
Eddington-limited. In the models including evaporation 
(see §1^) e varies to reflect the very low efficiency of the 
ADAF. Without evaporation, e is constant, equal to 0.1. 

C is a measure of the fraction of the X-ray luminos- 
ity that heats up the disc and, as such, contains infor- 
mation on the irradiation geometry. X-ray albedo. X-ray 
spectrum, etc. Kim et al. ( |1999 ) used an analogous pre- 
scription for their indirect irradiation flux. We found that 
the observed optical magnitudes and stability properties 
of persistent low mass X-ray binaries were compatible with 
a value of C« 5x10"^ (^j^i^ e = 0.1, Dubus et al. [1999[ ). 
Esin and co-workers found some evidence for changes of 
C in the outbursts of GRO J1655-40 (a.k.a. V1033 Sco), 



A0620-00 (a.k.a. V616 Mon) and GRS 1124-68 (a.k.a. GU 



ing of the outer disc in long period SXTs wou ld pr oduce 
a varying value of C (e.g. Ogilvie & Dubus, 2000). For 
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simplicity, we have nonetheless restricted ourselves to C 
constant in this study. 

2.4- Temperature dependence for a 

As in our previous works, we assume that a takes on two 
different values on the hot and cold branches of the 5- 
curve. Using a single value for a produces only low ampli- 
tude outbursts as numerous previous works have shown. 

Since irradiation modifies the hot branch, we adopt 
a slightly different prescription for the dependence of a 
on te mperature than that used by Hameury et al. ( 1998 ; 
1999|): 



log(a) = log(ac 



[log(ah) 



-log(ac)] 

rri \ 8 

^ crit 

~t7 



(9) 

where Tciit — 0.5 [T'c(Sniax) +20(2111111)] is calculated using 
Eq. I and |. This ensures that, even with irradiation, a — 
ah on the hot branch and a = etc on the cold branch. 

2.5. Boundary conditions 

The b oundary conditions are those of Hameury et al. 
( |l998|) . We stress that in all models, the outer disc ra- 
dius i?out is variable. The removal of angular momentum 
from the disc by the tidal fo rces fr om the secondary (as 
expressed in Hameury et al. 1998 ) sets i?out- The mass 
transfer rate Aftr from the secondary to the accretion disc 
appears in the outer boundary conditions. 

The inner boundary conditions change when evapo- 
ration is included and this is discussed below in §4. In 
models without evaporation, the inner disc radius is held 
fixed at a given i?in. In our calculations we chose larger 
values than appropriate if the disc were to extend down to 
the compact object e.g. to i?in = 3i?s for a black hole pri- 
mary. This is essentially a numerical limitation due to the 
very low densities at low radii when the disc enters quies- 
cence (see Fig. ^). This is not a drawback: as stressed in 
the introduction, there are strong arguments in favour of a 
truncated disc in quiescence rather than a disc extending 
down to the compact object. 

In outburst i?in could be expected to be close to the 
compact object but our assumption does not change any 
of the results. In outburst the inner disc is hot and the 
viscous timescale at which density perturbations evolve 
is very short at low radii. Adjustments of the inner disc 
are therefore instantaneous compared to the evolution 
timescale as long as i?in is low enough. We have verified 
this by computing some models with i?i„ = 10^ - 10^ cm 
which showed identical outburst lightcurves. A large i?in 
also avoids the Lightman-Eardley instability which arises 
when the radiation pressure dominates in the inner disc. 
Note however that str ong ir radiation can suppress t his in - 
stability (Czerny et al. 1986 ; Mineshige & Kusunose 1993|) . 



2.6. Numerical considerations 

The adaptive mesh numerically resolves the heat fronts 
and this is one of the major advantages of this scheme 
(Hameury et al. 1998| ). The disc mass is not necessarily 
conserved in the method and therefore provides a test of 
the numerical accuracy. At least 1200 mesh points were 
used in the runs shown here, resulting in excellent mass 
conservation for all cases. The only previous attempt at 
including irradiation heating in the DIM used at most 21 
radial zones, casting very serious doubts on the validity 
of their results (Mineshige et al. 1990; Kim et al. 1999) 



Hameury et al. ( 1998[ ) have emphasized the importance 
of having numerically resolved heat fronts (i.e. requiring 
much more than 21 grid points) to obtain reliable results. 
In another study, Cannizzo (1998) calculated the expected 
irradiation temperature throughout the eruption cycle of 
an unirradiated disc but not its feedback on the evolution. 

As starting points to our calculations we used the (un- 
stable) steady-state solutions. We then let the disc relax 
to the periodic outburst-quiescence cycle which is inde- 
pendent of the initial starting point (see §4.7 of Menou et 
al. [I999a|). 



3. Non-irradiated discs 

The evolution of the inner mass accretion rate times the 
efficiency gives a measure of the bolometric lightcurve of 
the system (Lx ~ eMnC^) (one should in principle con- 
volve this with the detector spectral response to compare 
with the measured X-ray light curve). In the following, 
we usually show only the evolution of Min, so these are 
not strictly speaking lightcurves, even though they have 
the same general characteristics. We show some typical 
'lightcurves' predicted by the standard DIM for parame- 
ters appropriate to soft X-ray transients before turning to 
models including irradiation. 

3.1. Re flares 

A characteristic feature of the DIM applied to SXTs are 
the multiple reflares that occur during the outburst de- 
cay (F ig. [1]) . These were previously discussed by Menou 
et al. (2000). Reflares appear when the surface density E 
behind the cooling front is high enough to reach Smax- At 
the radius at which this happens, the disc becomes ther- 
mally unstable and a new heating front develops. This 
front propagates outwards like an inside-out outburst, re- 
heating the disc until S(i?) ^ Emin, when cooling can re- 
sume. The density in the cold region is depleted as matter 
is accreted during this process, and the following reflare 
occurs at smaller radii and have lower amplitudes (Fig. ^ . 

Menou et al. (1999a) found the density behind the 
cooling front is roughly if Emin where K is a function of 
the v iscosity a and of the primary mass Mi (Vishniac 
7 M0 , if w 6 — 7 so that if Smin can be- 
and trigger a heating front. K is 



19971) . For Mi 



come greater than E„ 
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3.2. A deficiency of the DIM 



10 



10' 



16 



f- 



400 600 
days 




200 
days 



Fig. 1. Typical SXT variations of the inner accretion rate 
predicted by the standard DIM. This shows the mass ac- 
cretion at the inner edge of the disc around a 6 M0 black 
hole (top) and a 1.4 neutron star (bottom). The disc 
is not irradiated and the inner radius is kept fixed at 
i?in = 5 X 10^ cm (top) and 2.5 x 10* cm (bottom). 
The other parameters are those of the 'A0620-00' black 
hole a nd 'Aql X-1' neutron star models of Menou ct al. 
( I2OOOD . The standard DIM apphed to SXTs is character- 
ized by short recurrence times, low accreted mass during 
outburst and the presence of strong reflares during the 
decline. Clearly, this is does not reproduce observations. 



lower for lower primary masses {K « 4 for Mi = 1.2 M©) 
making it less likely to have reflares. Menou et al. noted 
reflares were absent in their calculations with a neutron 
star primary. 

Reflares also appear for neutron stars, but at smaller 
radii than considered by Menou et al. The reason is that 
the post-front S is not strictly proportional to Emin and 
has a slightly shallower dependence on radius: S oc R^-^ 
for a NS (Fig. |) instead oiR^-^ for and S^ax (Eq. |- 
^ without irradiation) . Hence there is always a radius for 
which S behind the cooling front will become greater than 
Smax and ignite a reflare. For a neutron star, the first re- 
flare appears after the cooling front has propagated well 
inside the disc (Menou et al. 2000, stopped their calcu- 



lations before this radius was reached) and has a strong 
effect on the lightcurve as it heats back a significant frac- 
tion of the disc (bottom panel of Fig. |^). 



Reflares are a generic feature of the DIM but do not cor- 
respond to observe d feat ures in the decay lightcurves of 
SXTs (Chen et al. 1997). Specifically, the rise-times are 
much too long in the calculated reflares. The decreas- 
ing amplitudes and timescales in-between outbursts is dif- 
ferent from the regular mini-outburst displayed by GRO 
J0422-h32 (a.k.a V518 Per). The number of reflares pre- 
dicted is also incompatible with the glitches of e.g. A0620- 
00. 

Reflares are therefore a deficiency of the DIM (Menou 
et al. 2000| ). Models of dwarf novae are not affected by 
this problem because the radii of the white dwarf primary 
(~ 10® cm) and of the discs (few times ~ 10^*^ cm) do 
not span a range large enough for reflares to appear. The 
reflares can be avoided in SXTs if, for instance, the thin 
disc does not extend much further in than the radius at 
which reflares ignite. However, a thin disc truncated at too 
high a radius will be cold and stable rather than transient. 

Menou et al. ( ^OOOD were able to avoid reflares by using 
a combination of evaporation (reducing their number) and 
low quiescent viscosity (increasing the density at which a 
heating front is ignited Smax)- Although the recurrence 
timescale predicted between outbursts was reasonable, the 
outburst shape did not compare well with FRED-type ob- 
served lightcurves, the rise time being too long. In addi- 
tion, the value for the quiescent viscosity was lower than 
usually assumed (ofc ~ 5 x 10"'^). This is not a problem 
per se considering the present knowledge on viscosity in 
quiescence (see e.g. Gammie & Menou 1998) but it could 
indicate different transport mechanisms in U Gem type 
dwarf novae with ac ~ 0.02 and SXTs with ac ^ 0.005. 
However, the mechanisms might be the same, but the ef- 
ficiencies different. However, very low values of ac have 



been advocated in the dwarf nova WZ Sge (Smak |1993D 
which is in many respects similar to SXTs; in this case, 
the very low value of ac must be due to a different mech- 
anism, since it is very unlikely that the efficiency of the 
same mechanism could vary by 2 or 3 orders of magni- 
tude. The following section will show how irradiation by 
lowering E behind the cooling front solves this problem 
without the need for very low values of a^. 



4. Irradiated discs 

In this section we study the effect of including self- 
irradiation to the standard DIM. We present only one 
model for this purpose (others may be found in Dubus 
pOOOl ) with Ml = 7 Mq, Mtr = IQi^ g s-\ ah = 0.2, 
ac = 0.02, < i?out >= 10^^ cm and a fixed inner ra- 



dius at i?in = 10 cm (this is further discussed below in 
§4.6). Irradiation is included using Eq. ^ with e = 0.1 
when Min > lO^^ g s"! and e = 0.1(Mi„/10i6 g s-'^f be- 
low. This amounts to a quick cutoff of irradiation below 
2q16 g s~^. The exact form of e matters little as long as it 
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Fig. 2. Reflares in non-irradiated discs around a NS (top 
panel) and a BH (bottom panel). The full lines arc the 
surface density E at different times during the decay from 
outburst (Fig. ^ and the dotted lines are Smax and Smin 
(Eq. H, H). In the top panel, the situation before the onset 
of the first reflare is shown by the top curve: S just behind 
the cooling front is close to the critical Emax for which 
the disc switches back to the hot state. The middle and 
bottom lines show the profile after several reflares: S in 
the outer disc has been depleted by the successive passages 
of the heat front (which can be seen in the bottom curve). 
The lower panel shows the situation after many successive 
reflares in a disc around a BH primary. The wiggles in 
the density profile corresponding to previous reflares have 
not yet been smoothed out by viscous diffusion. The inner 
edge is at 10^ cm (top) and 5 x 10* cm (bottom). 



ensures that irradiation is only important during the out- 
burst (a reasonable assumption). The resulting outburst 
time profile is shown in Fig. ^ 

Rise 

We start with a cold, quiescent disc in which most of the 
matter has been accreted in a preceding outburst. During 
quiescence (see below), mass transfer from the secondary 
replenishes the flow until the density becomes high enough 
and an annulus arrives at the thermally unstable branch 
of the S'-curve. The ring cannot maintain thermal equilib- 
rium and undergoes a transition to the hot branch. 

As in the standard DIM, two heating fronts are formed 
which propagate inwa rds and outwards from the ignition 
radius (Menou et al. |l999a and references therein). In 
'inside-out' (type B) outbursts the ignition radius is small 
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Fig. 3. Example of an outburst cycle when irradiation is 
included. From top to bottom: the mass accretion rate at 
the inner edge Min (full line) and Mi„ = eMin (dotted 
line), the V magnitude, the outer disc radius i?out (full 
line) and the transition radius -Rtians between the hot and 
the cold regions (dotted line), the mass of the disc Mdisc- 
The parameters are Mi = 7 Mq, Mtr = 10^^ g s"\ ah = 
0.2, Qfc = 0.02, < i?out >= 10^^ cm and a fixed inner 
radius at Rin — 10^ cm. Details of the outbursts and the 
evolution of the density and temperature profiles can be 
seen in Figs. ^[T 



and the propagation time of the inward bound front is 
very small compared to the other one. In 'outside-in' (type 
A) outbursts, the opposite is true. In the standard DIM, 
'inside-out' fronts can stall rather easily leading to short, 
low-amplitude outbursts. The density profile in quiescence 
is generally not far from Smax cx: i? so that, as the out- 
ward front progresses through the disc, it encounters re- 
gions of higher densities and Smax (the critical density 
needed to raise a ring to the hot branch). If the front does 
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not transport enough matter to raise the density at some 
radius above Smax, it stalls and a cooling front develops. 
'Outside-in' fronts always progress through regions of de- 
creasing Einax and therefore always heat the whole disc. 
Another consequence is that inside-out fronts propagate 
slowly, leading to slow rise times. This is particularly evi- 
dent in Fig. 8 of Menou et al. ( |200C| ). 

Irradiation does not change the structure of the heat- 
ing fr ont (compare Fig. ^ with Fig. 1 of Menou et al. 



aoo 



1999a) but does change the maximum radius to which 
an inside-out outburst can propagate. Since we do not 
consider self-scrccning, the outer disc always sees the ir- 
radiation flux from the central regions. As an inside-out 
front propagates, Mjn rises and the outer cold disc is in- 
creasingly irradiated (see Fig. |^). Irradiation heating re- 
duces the critical density needed to reach the hot branch 
(Eq. ^), easing front propagation. Obviously, a larger hot 
region implies a greater optical flux and irradiation always 
lowers the peak optical magnitude. 

In an inside-out outburst, the outer regions the front 
has not yet reached are frozen on the timescale on which 
Mill evolves (see Fig. |]). If Mjn (hence Mirr) increases on 
timescales shorter than the thermal timescale in the cold 
disc this can lead to situations where Ti„ > Tc at some 
radii. The vertical heat flux changes sign and would re- 
quire negative values of a in our treatment of thermal im- 
balance i.e. the assumptions of the code break down. We 
assumed that such an annulus is dominated by irradiation 
i.e. is isothermal at Tc and that irradiation contributes 
an additional heating term to the radial thermal equation 
(Eq. 0) Q+ ^ ^ cr(7]^, - r<f) to reflect the imbalance at 
the photosphere between the outgoing flux aT^ and the 
incoming flux crT^.^. This, or other assumptions, actually 
had very little influence on the rise-to-outburst lightcurve. 
Further studies of this phenomenon would require a de- 
tailed model of irradiation in SXTs where the geometry, 
the exact flux and spectrum of irradiating photons are 
properly set out. Note also that the absence of any back 
loop in the equations to prevent this situation may sug- 
gest that some additional physics is needed. One possibil- 
ity is that rapidly increasing irradiation woul d eva porate 
the u pper layers of th e disc (Begelman et al. 1983 ; Hoshi 
19841 ; Idan & Shaviv |l996| ; de Kool & Wickramasinghe 



199C) 



Fig. 1^ shows the evolution of the S and Tc radial pro- 
files during the outburst rise. As the front reaches the 
outer edge of the disc, the profiles in the hot region con- 
verge to those of a steady disc with constant M{R). This 
is because the viscous timescale, which is inversely pro- 
portional to aTc, becomes short enough to equilibrate 



50 



-50 



1 ' 1 ' 


T ' 1 ' 


10" 


1 ' 1 ' 1 ' 1 ' 








T.„ \ 






10^ 




s / 






{ 






\ 


— 1 — 1 — 1 — 1 — 1 — 1 — 1 — > 


1 1 1 1 1 1 1 


\ 




0.2 


a 


1 


0.1 




1 , 1 , 1 , \ > 








1.30 1.35 
R (10"' cm) 



1.30 1.35 1.40 
R (lO"" cm) 



Fig. 4. An 'inside-out' heating front in an irradiated disc. 
This is the front at i « 1.5 days in Fig. ||. S is in g cm~^, 
temperatures are in K, M is in units of 10^^ g s~^. The 
cold outer disc is almost isothermal with Tc « Tin.. 



disc regions does the vertical structure become irradiation- 
dominated, i.e. isothermal. 

The peak accretion rate (and optical magnitudes) is 
rather difficult to estimate analytically since it depends 
on the maximum radius to which the heat front can prop- 
agate. Qualitatively, this is set mostly by the conditions 
in the disc at the onset of the outburst. The principle fac- 
tor is the total disc mass when the outburst starts Mmax 
which is constrained by Smax- Parameter studies (§6) show 
that any changes which result in a lower Smax, hence in a 
lower Mdisc do lead to smaller outburst peaks. Increasing 
Ml, Qfc or decreasing i?out give smaller Mpcak- Increas- 
ing the mass transfer rate leads to higher disc masses in 
quiescence and higher peaks. 

In any case, Mpoak should be lower than or of order 
of the critical mass accretion rate Merit corresponding to 
Smin at i?out, as the disc is close to steady state at maxi- 
mum if it has been brought entirely to a hot state. If the 
heating front has not been able to reach the outer edge, 
Mpcak will be lower. 



the mass flow in the hot region (e.g. Menou et al. 1999a, 
and references therein). Irradiation has little influence on 
the actual vertic al str ucture in this region as discussed 
in Dubus et al. (|l999|) and we flnd Tc oc E oc i?-^/^ as 

in a non-irradiated steady disc. Only in the outermost Afpeak oc (ofh/Qc)*'"^. 



^ The dependence is really on the ratio ah/cvc. The mass 
in the disc is almost the same between the beginning and 
the peak of the outburst. Assuming Mmax oc Emax oc a^" ® 
and Mmax oc M^'J^^^a^"'^ (Shakura-Sunyaev steady disc) shows 
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4-. 2. Decay from outburst maximum 

The outburst decay can generally be divided into three 
parts: 

— First, outer disc X-ray irradiation inhibits the cooling- 
front propagation. Since more mass is lost than gained 
from the secondary {Min 3> A/tr) the disc is drained 
by viscous accretion of matter (King & Ritter 1998). 

— Second, the accretion rate becomes low enough that 
the X-ray irradiation is unable to prevent the cooling 
from propagating. The propagation speed of this front 
is, however, controlled by irradiation. 

— Irradiation plays no role and the cooling front switches 
off the outburst on a local thermal time. 

We now discuss these three phases in more detail. 
4.2.1. Viscous "exponential" decay 

In the model presented here the disc becomes fully ionized 
at the outburst peak. The disc then evolves with M{R) 
almost constant so that vYj ^ M\a{t) /"in] the total mass 
in the disc is thus 

Afd = j 2-KRY.dR oc A/in j \-dR (10) 

and one can consider that the disc decays through a 
sequence of quasi-stationary states. In irradiation domi- 
nated discs as well as in standard Shakura-Sunyaev discs, 
u(xT (X MP/^^+P)^ with (3 = 3/7 in hot Shakura-Sunyaev 
discs, and /3 — 1/3 in irradiation dominated discs. The 
outer disc radius does not vary much so the time evolu- 
tion of the disc mass is: 



dA/d 
"dT 



= -A/in OC M 



1+/3 



(11) 



showing that A/in evolves almost exponentially, as long 
as A/^ can be considered as constant (i.e. over about a 



decade in A/i„, as found by King fc Ritt er 1998 for irradi- 
ated discs, and by Mineshige et al. 1993| for non-irradiated 
discs with angular momentum removal). Power-laws with 
indices close to -1 are found in discs with constant angu- 
lar momentum, i.e. in which there are no tidal torques 
preventing the outer di sc ra dius to expand i ndefin itely 
(Lyubarskii & Shakura |l987l Cannizzo et al. |1990|; Mi- 



neshige et al. 1993), and are of no practical interest here. 

The decay is viscous as long as thermal equilibrium can 
be maintained. Incidentally, this makes it difficult, if not 
impossible, to have a viscous decay in the standard DIM 
without additional assumptions. In outside-in outbursts, 
matter accumulates at the outer disc edge and the result- 
ing density excess decreases viscously, producing 'flat-top' 
light curves. Thermal equilibrium requires the whole disc 
to be kept in the hot state. If the outer disc is cold then 
there will be a thermally unstable region evolving rapidly 
(i.e. non- viscously, see next paragraph). Keeping a disc 
hot can be achieved by additional mass transfer during 
the outburst or non-standard sources of heating (e.g. tidal 



heating as in Buat-Menard et al. 2000| ). A smaller ac in- 
creases the difference between Emax and Smin and hence 
can also help keep the disc in the hot state longer (Menou 
et al. 2000). 

Irradiation heating not only leads to the exponential 
decay but also provides a natural way to keep the disc 
hot. An irradiation temperature at the outer edge above 
10'* K ensures hydrogen is everywhere ionized. For such 
temperatures, the disc is thermally stable whatever the 
local density S and accretion rate. Put differently, the 
lower branch of the iS-curve disappears when Tirr > 10"* K 
(see Fig. 4-5 of Dubus et al. 1999) or S,„in — Emax at 
Tirr w 10'* K (and are undefined below: see Eq. ||, ||). 

One should remember, however, that many SXT light- 
curves are not exponential despite their discs being irra- 
diated. Clearly additional physical processes have to be 
taken into account if one wishes to explain such 'non- 
typical' behaviour (see e.g. Esin et al. ^OOOa ). 



4. 3. Irradiation- controlled linear decay 

The second stage of the decay starts when a disc ring can- 
not maintain thermal equilibrium and switches to the cool 
lower branch of the S'-curve. The disc outer edge always 
becomes unstable first with E < Emin since E oc i?~^/'* 
(see above) while Emin oc R. In an irradiated disc this 
happens when the central object does not produce enough 
X-ray flux to keep the T'irr(/?out) above 10'* K so that 
the low state is again possible. The material in the low 
state has a lower viscosity and piles up, leading to the 
appearance of a cooling front with a width w cc H (Pa- 
paloiz ou fc Pringle |1985| ; Meyer |1986| ; Fig. 4 of Menou 
et al. 1999a). The front propagates through the disc at a 
speed Vfiont ^ OhCs (Meyer 1984 ; Vishniac & Wheeler 



1996). In the standard DIM, the sound speed depends on 



the temperature at the transition between the hot and 
cold regions and is almost con stant a s verified by numeri- 
cal calculations (Menou et al. 1999a ). 

In an irradiated disc, however, the transition between 
the hot and cold regions is set by Ti^r since the cold branch 
only exists for Tirr < 10'* K. Unlike a non-irradiated disc, 
the cooling front can propagate only as far as the radius 
at which Ti„ w 10'* K, i.e. as far as there is a cold branch 
to fall onto. Here also irradiation controls the decay. The 
hot region stays close to steady-state but with a shrinking 

size /?hot ~ (Eq. | with Tir,(/?hot) - const). The 

mass redistribution during front propagation is complex 
with the hot region losing mass both through inflow and 
outflow and our models show steeper declines than what 
some analytic approximations predict (e.g. King 1998). 

The linear decay is illustrated in Fig. |^. The irradia- 
tion temperature is lower than 10* K at the outer edge for 
t > 70 days where the decay becomes steeper. The cooling 
front (see also Fig. ||) appears as a depression propagating 
inwards in the surface density profile. The irradiation tem- 
perature and the critical densities Emin.max are plotted for 
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Fig. 8. The cooling front in an irradiated disc. This is the 
front at t « 190 days in Fig. |[ E is in g cm~^, temper- 
atures are in K, M is in units of 10"'^^ g s~^. The front is 
at the position for which Tirr (dotted Hne) is « 10^ K. At 
this point Smin « S,nax sincc there is no cold branch for 
higher Tjir (smaller R) . As Tirr decreases the two critical S 
separate and converge to their non-irradiated values (this 
can be seen in Fig. |^; note the values of Smax and Emin 
shown here are only accurate to the precision of the fits 
given by Eq. |,|). The cold outer disc is almost isothermal 
with Tc w Tii-i-. The negative values of M show the cooling 
front transports matter to the outer disc. 



the last profile at t « 190 days, showing the cooling front 
is at the radius for which Tiir w 10^ K and that Smin,max 
are undefined for lower radii. 



4.4- Final thermal decay 

The quick decay after t > 190 days is due to vanishing 
irradiation as e becomes very small for Afjn < 10^^ g s^^. 
The cooling front thereafter propagates freely inwards, on 
a thermal time scale. 



4-5. Mini-reflares 

As discussed above, the edge of the hot zone is at the 
stability limit set by Tirr ~ lO** K. At this point, one has 
E « Emin « Einax bccausc irradiation modifies the critical 
surface densities. Further away from the edge of the hot 
zone, the irradiation temperature decreases rapidly {T-^^ ck 
i?^^) and Eniax gradually becomes greater than E„iin. In 
the non-irradiated limit, the ratio Emax/Smin is a constant 
depending only on the ratio ac/a^ (see e.g. Fig. ^). 




Fig. 9. Example of mini-reflares in an irradiated disc. 
The region close to Tirr ~ 10^ K is highly unstable with 
Emin w Emax (Fig- H) • This Can lead to a succession of 
heat fronts propagating quickly back and forth in a small 
region (bottom panel). The exact details depend on the 
model assumptions (e.g. a). This has no influence on ei- 
ther Min (top panel) or the optical flux (not shown). See 



f;4.5 for details. 



The region immediately behind the cooling front where 
Smax is close to Emin is clearly very unstable. Slight vari- 
ations of E may suffice to have E > Emax behind the 
cooling front and therefore start a reflare just as in a non- 
irradiated disc. This depends on the numerical details of 
the model, such as the functional of a with temperature 
(§0). 

These mini-reflares are unimportant. The critical den- 
sity Emin increases quickly behind the cooling front (see 
Fig. ^) because the irradiation flux decreases with 
Therefore, the heating front in a mini-reflare reaches al- 
most immediately a radius where E < Emin and cooling 
resumes. Furthermore, irradiation seriously depletes the 
disc and the front flnds very little matter to fuel its prop- 
agation. Ejiiiii is much lower than its non-irradiated value 
(Fig. ^). For instance, Fig. 4 of Dubus et al. ( 1999| ) shows 

Smin(Tirr « lO^K) « 30 g Cm'^ instead of E,rrin(Tirr = 

OK) « 150 g cm~^. The much lower post-cooling front 
densities prevent large reflares as in the standard DIM. 
In practice, mini-reflares in irradiated discs have no influ- 
ence on the lightcurve (very low amplitudes and cycles, 
see Fig. |) but can be a numerical nuisance. 

4-6. Quiescence 

As discussed above, irradiation dramatically slows down 
the cooling process. The decay lasts much longer and, ac- 
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cordingly, the total mass accreted during an outburst is 
much larger. The densities after propagation of the cool- 
ing front are low so the total mass of the disc as it enters 
quiescence is much less than in the non-irradiated case. 

The low quiescent surface densities imply low midplane 
temperatures. The (5],Tc) S'-curves show that for low S 
and no irradiation (as appropriate once the disc enters 
quiescence) the lower branch is flat with Tc « 2000 K 
whatever R or Mi (see e.g. in Figs. 4-5 of Dubus et al. 
1999). In contrast, the standard DIM leads to higher S, 
close to the non-irradiated Emin and thus to Tc w 4000 K 



(e.g. Osaki |1996| , Menou et al. ^OOOD . 

The critical density Emax that the disc must reach be- 
fore an outburst can start will be the same as in the stan- 
dard DIM if irradiation in quiescence is negligible. This is 
in general true since Min oc X^^. is very smalP. Since the 
quiescent disc after an irradiation-controlled decay is al- 
ways less massive, the time to build up the critical density 
will be longer. 

In the model presented in this section. Tin- = in qui- 
escence. The evolution of the E and Tc profiles (Fig. ^ 
shows matter arriving from the secondary at the outer 
edge diffuses down the cold disc. The disc contracts due 
to addition of this lower angular-momentum material. The 
temperature is almost constant in most of the disc with 
Tc « 2000 K. This is not true in the later stages of qui- 
escence where Tc rises in the inner disc due to the build 
up of matter there: the inner rings move up along the 
lower ranch of the (E,Tc) 5-curve and out of the flat 
Tc sa 2000 K region discussed above. In the last profile, E 
is close to Eniax (dotted line) and the outburst is ignited 
at i? « 2 X 10^ cm. 

The ignition radius depends essentially on the mass 
trans fer rate from the secondary and the disc's size (Smak 



198J). For low Mtr, matter drifts down the cold disc and 
Emax is reached in the inner disc (inside-out, type B out- 
burst). As Mtr increases, the mass accumulation time at 
the outer radius can become lower than the drift time, 
triggering an outburst in the outer disc (an outside-in, 
type A outburst). The accumulation time is (Ichikawa & 
Osaki |1994| ; see also Osaki |l995| , |l996| and Lasota |200lD : 

= 2^i^outEn,ax(i?out) ^ ^^2) 

Mtr 

A ring of material spreads viscously over AR « (j^tacc)^^^ 
during tacc- Setting iacc = ^drift gives an estimate of the 
accretion rate for which the transition from type B to A 
outbursts is expected. Using Eq. || and tdrift = St^i s {S is 
a numerical correction factor introduced by Osaki 1996 ) 
outside-in outbursts are therefore expected when 
Mtr > 27ri?„„t5-i/2i/E, 



> 3.3 X 10 



16 



a 
002 



(^out) 
"■2 / Ro 



10 



11 



cm 



2.6 



^ The qui esce nt disc can be irradiation-dominated if C is very 
In this case, the critical Emax is lowered by 



large (see §6.1 



Ml 



-0.9 



f Tc 



7MqJ V2000Ky ^ ^ 
with 6 = 1. This is smaller than the mass accretion rate 
needed to stabilize an irradiated disc with these parame- 
ters (Mtr Z 10^^ g s-\ Eq. 30 of Dubus et al. |1999D. In 



principle type A outbursts are thus possible in SXTs, but 
numerical calculations show that real outside-in outbursts 
starting far out in a large disc are difficult to obtain. Osaki 
( 1996D finds S w 0.05 provides an adequate fit to his cal- 
culations. The quiescence time of our model (about 4000 
days) implies idrift ~ O.ltvis(^out) so that the required 
A/tr for a type A outburst is above the accretion rate for 
which the disc is stable. 

The inner edge of the accretion disc plays a crucial 
role in quiescence. For low accretion rates Emax is reached 



needed to trigger an outburst is reduced when the disc 
inner radius is smaller: the higher i?in in quiescence, the 
longer the recurrence time. Even with = 10^ cm the 
model is still short of providing 

— the long recurrence timescale (ircc ~ 10 years com- 
pared to the tens of years inferred in SXTs) and 

— the quiescent X-ray luminosities: the predicted quies- 
cent accretion rate, in the range 

1012-13 g g-i is gyn 
too small to account for observed X-ray luminosities 
exceeding 10^^ erg s~i. 

The next section will show that including disc evapora- 
tion into a hot, low-density, accretion fiow can solve these 
problems. 



5. Irradiated discs with evaporation 



Following the work of Menou et al. ( 2000 ), we now assume 
that the inner disc is gradually evaporated into a hot, 
radiatively inefficient accretion flow during the decay from 
outburst. The exact nature of the flow is not important 
as long as it does not participate in the dynamics of the 
thin accretion disc. A good assumption is that the inner 
disc is replaced by an advection-dominated accretion flow 
(ADAF). ADAF-fthin disc models have been successful in 
explaining the spectral states of SXTs and in particular, 
the quiescent X-ray flux provided the mass accretion rate 



in the ADAF is large enough (e.g. Esin et al. 1997). 

The accretion timescale in an ADAF becomes of the 
order of the thermal timescale so that the inner flow can 
be considered as quasi-steady and does not participate in 
the limit cycle. A good assumption in our dynamical study 
is therefore to treat evaporation as a variation of the inner 
radius of the thin disc we model. We take the same rate 



of evaporation as in Menou et al. (2000) 



Mcv(i?) = 0.08A/Edd 



1/4 

^1 +£ 



irradiation and the recurrence timescale is shortened. 



VSOOi?, 

The inner radius of the disc is defined as: 

Mev(-Ri„) = Mi„ 



n -1 



(14) 



(15) 
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Below this radius all the matter evaporates. Because of 
the steep dependence of Mev on R, we can safely neglect 
evaporation above this radius. As noted in Menou et al. 
( 2OOOD , the detailed functional dependence of Mev(-R), or 
of i?in(Min) has essentially no effect on the results; what 
matters is the value of the inner disc radius during quies- 
cence. 

Equation is entirely ad hoc, and is not based on 
any particular physical mechanism. There are models that, 
in principle, allow for the determination of the evaporation 
rate as a function of radius. However, the variety of mod- 
els shows that there is no general agreement on the phys- 
ical cause of evaporation. For example, Meyer & Meyer- 
Hofmeister (1994) assume that electron conduction from 
the disc to the corona plays a major role, whereas Shaviv, 
Wickramasinghe & Wehrse (1999) consider a thermal in- 
stability related to the opacity law. In addition, even for a 
given model, the evaporation law depends on parameters 
which are not easily measured (as for example the mag- 
netic field and its coherence in the case of electron con- 
duction ; Meyer, Liu & Meyer-Hofmeister 2000). As all 
these models have been proposed to explain the existence 
of holes in observed accretion discs, they, by definition, 
tend to produce similar results, i.e. the truncation radius 
must be at a detectable distance from the compact objet. 
And of course, for the same reason, Eq. ( |l^ would also re- 
produce these results. The latest calculations by Meyer et 
al. (2000) give evaporation rates of the order of 1.4 xlO^^ 
g s~^ at 10^°cm from a 6 Mq black hole, and 1.6 xlO^'' 
g s~^ at 10^ cm. These are larger by a factor ~ 15 and 
~ 2 respectively than the values we are using here. Note 
that the evaporation rate is so high in Meyer et al. (2000) 
that the quiescent disc would be truncated at too large a 
radius to be unstable. 

The ADAF efficiency roughly sc ales a s M (i.e. the lu- 
minosity scales as M^, Esin et al. 1997 ) . From Eq.|l^, 
this is equivalent to e oc R^- In the model, we as- 
sume e = 0.1 when i?in = i?,nin (in outburst) and e = 
0.1 X (i?min/^in)~^ whcu i?in > -Rmin (in quiesccncc). Ir- 
radiation is usually negligible in quiescence so these as- 
sumptions have little importance there. Varying e may 
change slightly the end of the outburst when evaporation 
becomes important. 

We took i?min = 5 X lO'* cm in all the following models. 
As discussed below in §^]^, the only time Rmin is reached is 
in outburst where the inner disc is in thermal equilibrium 
and close to steady-state. The inner disc edge plays no 
role in this case. 

In the following we discuss in some detail the same 
model as in §4 but including evaporation. The overall 
outburst time profile is shown in Fig. |l^. A compari- 
son with Fig. ^ shows the outburst is almost identical 
but that the disc spends much more time in quiescence 
(tree ~ 21 years). The maximum luminosity reached dur- 
ing outburst is also increased. 



10 



IB 



17 



10 



io'5 

IS 

20 

> 25 

30 

35 
10.0 

^ 1.0 
5 

0.1 
10 

B 

DO 

2 ^ 
4 















rs — 




1 1 ! 




" ■ 









10 



so 

years 



30 



40 



Fig. 10. Example of an outburst cycle when irradiation 
and evaporation are included. The parameters of the 
model are the same as in Fig. |[ From top to bottom: Mjn 
(full line) and Afin- (dotted line); V magnitude; i?out (fuU 
line), i?trans (dotted line) and R^ (dashed line), Afdisc- 
Details of the outburst and the evolution of the density 
and temperature profiles can be seen in Figs. |ll[|l^. 



5.1. The outburst 

After an outburst is triggered, the mass accretion rate 
at the inner edge gradually increases as in the previous 
model. The inner radius of the thin disc decreases until 
its minimum value is reached for = Afcv(fimin)- A 
slight change of slope is associated with this (at the 4th 
dot in Fig. |ll|). The time t^-isc the disc needs to reach Rmin 
is a viscous timescale. Before the onset of the outburst, 
i?in w 6 X 10^ cm. Therefore, tj-iso is of the order of the 
difference between fvis at 6 x 10^ cm and 5 x 10* cm. The 
viscous timescale is given by 

Ui. = ^ = (GM,i?)V^^. (16) 
u akin 
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For parameters appropriate to this case {Tc « 10^ K, 
a = 0.2, fi = 0.5, Ml = 7 Mq) the difference gives about 6 
days which is in good agreement with the numerical cal- 
culation. Fig. |l^ shows the minimum radius is reached at 
t « 5.5 days, about 3 days after the onset of the outburst. 

The subsequent evolution is identical to the model 
where the inner radius is kept fixed. The peak accretion 
rate and optical flux are higher in the model with evap- 
oration because of the different quiescent histories of the 
discs. When the outburst starts, the discs do not have the 
same density distribution and total mass (Fig. |[ |lO|). 

The decay from outburst is also identical to the previ- 
ous model with a viscous and a linear decay (Fig. |l^) . The 
linear decay has a limited influence on the lightcurve since 
the cooling front only propagates to 8 x 10^° cm before 
evaporation sets in (t « 160 days in Fig. ^2|). Irradiation 
then shuts off with the decreasing e and the front cools the 
whole disc quickly. At the same time, the inner disc radius 
rises with the increasing evaporation which also quickens 
the cooling. In general, evaporation and/or changes in the 
irradiation efhciency e take place for Mjn ^ Mtr at which 
point the total disc mass is close to its minimum value 
(the disc starts filling in again) and there is little mat- 
ter left to accrete. Therefore, the exact prescriptions for 
.^evap and e do not have a significant impact on the out- 
burst lightcurve. For instance, these produce a difference 
of only a few days in the outburst lengths between the 
model of M and the model of this section. 



5.2. Quiescence with evaporation 

The whole disc becomes cold and enters quiescence when 
the cooling front reaches the varying inner radius Rin. 
In Fig. Il^ this happens at t w 195 days when i?in « 
5 X lO^^cm and Mi„ « 5 x lO^^ g g-i_ xhe drift time 
of cold matter is long and the disc cannot maintain this 
high accretion rate: Afin decreases as material is gradually 
accreted. This adjustment happens on a long timescale 
(cold material) which depends on the S profile left after 
the outburst. With a fixed i?in, Min is already very low 

(about 10^2 g g-i^ Pig |) .^yi^en the whole disc becomes 6. Parameter study 



very large disc sizes of SXTs, the time it takes for matter 
to diffuse down the disc is shorter than the accumulation 
time at the outer disc. A truncated disc does not change 
this conclusion: evaporation does not prevent inside-out 
type B outbursts. 

However, evaporation does suppress the small mini- 
outbursts which are found when the inner disc radius is 
small (see e.g. Hameury et al. 1998). With i?in fixed at 
10^ cm, Smax is small, of the order of 10 g cm^^, and 
an outburst can be triggered easily as soon as matter dif- 
fuses in to the inner radii. This leads to sequences of mini- 
outbursts which are not observed. A truncated disc will 
need to build up more mass before an outburst can start. 
The critical density Emax varies with R and will be much 
higher when i?in is large, of the order of a 100 g cm~^ at 
10^" cm. This prevents small outbursts from being trig- 
gered. 

A truncated disc at the onset of an outburst will be 
more massive than a non-truncated disc (compare Figs. |^ 
and H^). But since the disc in outburst cools under the 
same conditions (dictated by Ti„), the mass of the disc at 
the end of the outburst will be roughly the same in both 
cases. If AM is the mass accreted during the outburst 
then we have AMtrunc > AM„o trunc- 

The quiescence time is the time it takes to replenish 
the mass lost during the outburst: 
AM 

iquicsc — ~v 7v (1 ') 

Irradiation depletes the disc during the decay, increases 
AM and leads to longer tquicsc- In contrast to models 
in which the disc extends to low radii, truncated discs 
can have larger Min in quiescence implying longer tquiesc- 
Eq. |l^ also shows that the quiescence time depends on 
the mass transfer rate even for inside-out outbursts when 
Min in quiescence is a significant fraction of Mtr. In a disc 
with low quiescent Min, iquicsc for inside-out outbursts is 
of the order of idiff which is independent of Mu- (Osaki 
1996| ; Smak [199^ ). 



cold. This is actually a necessary condition for the disc at 
low radii to reach the cold branch. For such low accretion, 
rates the supply of material is enough for Afin to increase 
steadily throughout quiescence. 

This first stage lasts until material from the outer edge 
had enough time to drift to the inner edge and increase 
Min- In Fig. |l| this happens at t « 3000 days. The first 
three E radial profiles of Fig. |l^ show clearly mass from 
the outer edge diffusing inwards and reaching the inner 
edge at t fa 2500 days. The disc then steadily builds up 
mass until the combination of a decreasing i?in (hence 
Emax) and increasing E triggers an outbu rst (se e also Can- 
nizzo |l998b| , Meyer & Meyer-Hofmeister |1999D . 

Under most circumstances (§4^) the outburst will be 
of the inside-out type. The reason is that, because of the 



6.1. Changing the irradiation strength C 

We show in Fig. |lj the effects of a stronger or weaker ir- 
radiation on the outbursts. The strength of irradiation 
is set by the value of C and we show models for C = 
10~^, 5 X 10"'^, 10~^, 10~^. As irradiation grows stronger, 
Emin becomes lower (see Eq. ||) and the first effect is to in- 
crease the length of the outburst: more mass gets accreted 
before Emin is reached. Thus, the relative amount of mat- 
ter accreted during outburst, AAZ/Af, increases with C. 
The disc spends more time accreting in quiescence before 
it can build up enough density somewhere to reach Emax- 
Like E„iin, Emax also decreases with increasing irradi- 
ation (see EqJI). Therefore, one would also expect Emax 
to be reached more easily, hence to have shorter quiescent 
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Fig. 14. Effect of changing the irradiation par ameter C on 
the outburst time evolution and cycle (see j |6.l| for details). 
C/10-3 =0 (fuU line), 1 (dash-dotted), 5 (dashed), 100 
(dotted) . All other parameters are identical to those of the 
model described in §4. The non-irradiated curve shows no 
reflares in this particular case. This is in general not the 
case, even with evaporation: the same model with Afi = 
6Mo shows them. 



times. Usually the irradiation flux is low at the end of qui- 
escence so only large values of C can heat the disc enough 
to change Smax (see Fig. |l|). 

There are two competing efl^ects: (i) the increased 
amount of matter accreted during an outburst lengthens 
the quiescence time; (ii) unusually strong irradiation in 
quiescence lowers the critical density needed to reach an 
outburst, which reduces the quiescence time. This happens 
when C is large, in which case irradiation is no longer neg- 
ligible in quiescence (in contrast to the models of §4-5). 
The numerical models show that there is a Cmax for which 
the recurrence time is greatest, C,„ax ~ 5 x 10^^, i.e. the 
'standard' value for which persistent low-mass X-ray bi- 
naries are stabilized (Dubus et al. 199£ ). This is probably 
no more than a coincidence as Cmax certainly changes as 
a function of Mi, i?out, etc. 

The average mass in the disc <Af(j> decreases with 
stronger irradiation. <Af(j> depends mostly on the mass 
of the disc in quiescence. Since S < Smax in quiescence 
and Sniax decreases with irradiation, <Mi> decreases as 
well. Another consequence of stronger irradiation are the 
higher optical fluxes in outburst due to the higher outer 
disc temperatures. 




R (10'°om) 



Fig. 15. Density profiles S(i?) in quiescent evaporated 
discs before the onset of an outburst. The dotted line is 
Smax- For C = 5 X 10^'' irradiation is negligible in qui- 
escence and the critical Smax is the same as for a non- 
irradiated disc. For larger values of C irradiation becomes 
important, lowering Smax- The critical density the disc 
must reach before an outburst is significantly less and the 
recurrence time between outbursts is decreased. See §3.1 
for details. 



6.2. Changing the viscosity a 

The effect of changing ofh and ttc is straightforward. The 
a parameter sets the viscous time of the accretion flow 
(Eq. |l^). A disc with a low viscosity in quiescence dif- 
fuses mass slowly to the inner edge. Therefore, the time 
to reach Emax is longer when ac is lower (see top panel 
of Fig. [l^ ). This also modifies the peak accretion rate of 
the outburst. Similarly, a disc with a higher viscosity in 
outburst accretes mass more quickly and the decay time is 
shorter (Eqs. [lo|-pT|). The conditions for the disc to enter 
quiescence and the subsequent evolution are independent 
of ofh- Therefore, the models shown in the top panel of 
Fig. |l^ which have different ah but the same ac have the 
same recurrence time. 



6.3. Changing the mass transfer rate Mtr 

Next, we explore the effect of changing the mass trans- 
fer rate Mtr, keeping all other parameters identical. For 
inside-out outbursts, the time spent in quiescence with 
fixed i?in is the diffusion time of matter from the outer to 
the inner edge, which does not depend on A/tr- However, 
models including evaporation show a dependence on Mtr 
when Min in quiescence is not negligible compared to Mtr 
(Eq. |l^ . Fig. ^ shows the quiescence time shortens with 
increasing Mtr- 

Higher Mtr result in longer outbursts since (i) the com- 
panion provides more material during the outburst and 
(ii) the total disc mass at the onset of the outburst in- 
creases with Mtr. The conditions in the disc are the same 
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Fig. 16. Effect of changing the viscosity parameters on 
the outburst time profiles. Bottom: changing ac modifies 
the recurrence time (full line ac = 0.01, dashed line ac — 
0.02, both with ah = 0.1). Top: changing ah modifies the 
outburst time evolution (full line cth ~ 0.2, dashed line 
ah = 0.4, both with ac — 0. 02. The recurrence time is the 
same iroc ~ 20 years). See 



.2 for details. 



at the end of the outburst, regardless of the mass transfer 
rate (except, of course, at the outer edge) because cooling 
happens in the same way for discs of the same size. A ring 
at a given radius can only cool when Ti„ reaches 10* K 
and the surface density behind the front is roughly the 
same in all cases. Thus, even with different Mtr, the discs 
have almost the same total mass as they enter quiescence 
(Fig. 0). 

In a truncated disc the drift of material in quiescence 
will change with Aftr because the distance it has to travel 
is rather short. Therefore, the conditions at the beginning 
of the outburst depend on the mass-transfer rate. In par- 
ticular, the ignition radius increases with Mtr as could be 
expected: a larger amount of mass from the outer edge 
diffuses down the disc and can trigger the outburst at a 
larger radius (but as mentioned earlier outbursts starting 
at the outer disc edge do not occur). This also increases 
the total amount of mass in the disc at the ons et of the 
outburst and hence the peak accretion rate (see §4.1). 



For the high mass transfer rates the disc stays fully 
ionized for a significant fraction of the outburst, leading to 
a viscous decay. Our model with the highest Mtr show this 
is accompanied by significant variations of the outer radius 
which produce deviations from the expected exponential 
(lower panel of Fig. |l^) . 
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Fig. 17. Effect of changing Aftr on the outburst cycles. 
Top two panels : Min and Mdisc for Mtr/10^^ g s~^ =, 
0.5 (fuU line), 5 (dashed fine) and 50 (dotted fine). The 
bottom panels zoom on the outburst for the same mod- 
els, showing A/in (full line) and Mi„ (dotted line), i?out: 
i?in and the transition radius between the hot and cold 
regions (dashed line) . The other parameters of the models 
are Mi = 7 Mq, i?out « 2.5 x 10" cm, ah = 0.2, ac = 0.02 
and C = 5 X IQ-^. See jfj for details. 



6.4- Changing the disc size -Rout 

The disc size in the model depends on the strength of 
the tidal truncation term in the angular momentum con- 
servation equation and on the circularization radius. The 
ratios i?ciic/a and i?out/a depend only on the mass ra- 
tio q — M2/M1. Assuming q « 0.1, which is reasonable 
for black hole SXTs, we derive -Rciic/^out ~ 0.5 (Pa- 



paloizou & Pringle 1977). In the following, we changed 
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Fig. 18. Effect of changing i?out on the outburst cycles. 
Top two panels : Afjn and Mdisc for <i?out>/10^° cm — 
9.3 (dotted line), 25.7 (dashed line), 107 (full line). The 
bottom panels zoom on the outburst for the same mod- 
els, showing (full line) and Afi^r (dotted line), i?outj 
and the transition radius between the hot and cold 
regions (dashed line) . The other parameters of the models 
are Mi 7 Mg, M = 5 x lO^^ g s-\ ah = 0.2, = 0.02 
and C = 5 X lO'^. See Slel for details. 



the mean outer disc radius, keeping the above ratio con- 
stant. This is equivalent to simulating binaries with the 
same components but increasing orbital period Porb ~ 
4.6 (i?out/10" cm)3/2M-i/2 hours {q ^ 0.1). 

Increasing the disc size makes it more difficult for the 
heat front to reach the outer disc radius and fully ionize 
the disc. In small discs the decay is entirely viscous while 
larger discs show only irradiation-controlled linear decays 
(Fig. |l^). In the case of a viscous exponential decay the 



decay rate increases with smaller disc radius in accordance 
with Eqs. ^0[|ll|. 

The ignition radius of the outburst increases for 
smaller disc sizes since the diffusion timescale is shorter. 
Smaller discs are also less massive and, therefore, take less 
time to replenish than large discs for the same mass trans- 
fer rate: the time spent in quiescence is shorter. This also 



leads to smaller outburst peaks (§4.1) 



6.5. Changing Mi 

Finally, we vary the mass of the accreting object Mi. Ir- 
radiation is implicitely assumed to be the same since the 
prescription for Tj„ (Eq. ^ does not depend on Mi . Smax 
and Smin have some dependence on Afi (Eq. |3| and ||) but 
their influence on the outburst is negligible compared to 
the changes induced by evaporation. The strength of evap- 
oration depends strongly on Mi through MEdd and Rs in 
Eq. ||: Mcv oc Ml (this cannot be seen in Menou et al. 
2000 where by error the evaporation formula for neutron- 
star SXT outbursts uses Mi = 6 Mg). In outburst, this 
modifies the accretion rate at which evaporation sets in 
which in turn reduces the length of the outburst for lower 
Ml (see also Meyer & Meyer-Hofmeister 2000 who stud- 
ied extensively the Mi dependence of the instability cycle 
with evaporation, but without irradiation, and reached 
conclusions similar to those of Menou et al. 1999). 

In quiescence, the disc around a lower mass compact 
object is truncated at smaller radii which reduces the in- 
ner accretion rate M^^. Mass is accumulated more eas- 
ily (Eq. |l^), leading to shorter recurrence time (Fig. |l9|). 
The lower accretion rates in the evaporated disc in qui- 
escence do not necessarily translate into lower quiescent 
luminosities. In our assumptions, the increase in accretion 
efficiency implied by a lower i?in (e c>c R^^^) compensates 
for the lower Min. Plotting Afin- — eM^ instead of only 
Afin as in Fig. ^, we find M,,, w lO^""" g s'^ for both 
Ml = 1.4M0 and 7M0. The actual efficiency of a 1.4 Af© 
neutron star may be much higher in quiescence than what 
we assume here. It might actually be more correct to as- 
sume e is always 0.1 for the neutron star since the accreted 
matter will radiate at the surface. In this case, the lumi- 
nosity in quiescence deduced from Ai^in will be higher than 
for a black hole primary. However, a propeller effect might 
significantly reduce the actual mass accretion rate onto a 
neutron star in quiescence (Menou et al. 1999c) making 
the distinction between black holes and neutron stars less 
obvious. 



7. Discussion 

Our most basic hypothesis is that the outbursts of SXTs 
are due to the thermal-viscous instability in a thin accre- 
tion disc. There is little doubt that LMXBs do possess thin 
accretion discs at least far from the compact object, while 
the wide applicability of the thermal-viscous instability in 
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suited not so much from the assumptions on the viscos- 
ity but on the physics which was or was not included. 
In CVs, studies that successively included a variable disc 
size, evaporation of the central regions, irradiation by the 
white dwarf, stream-impact and tidal dissipation heat- 
ing have gradually led to increased improvement in the 
agreement between observations and detailed predictions 
of models (the most recent being Hameury et al. |1998; 



1999| ; |200q ; Buat-Menard et al. |200(]| ). For example, it was 



Fig. 19. Effect of changing Mi on the outburst cycles of 
the inner mass accretion rate Mi^. The model of §5 is 
shown by the full line (Mi = 7 Mq). The same model but 
with All = 1.4 Mq is shown by the dotted hne. See § |6.5| 
for details. 



CVs gives strong support to its existence in LMXBs. It 
has long been known that models of dwarf novae require 
a lower viscosity in quiescence than in outburst. Models 
with ttc = ah predict low-amplitude outbursts with short 
cycles w hich d o not agree with observations (Smak 1984 ; 
Lin et al. 1985 ). The upper hot branch of the local thermal 
equilibrium curve (5'-curve) is characterized by a high vis- 
cosity while the opposite is true for the lower cold branch. 
Our 'standard' model refers to the thermal- viscous insta- 
bility in such a two-a thin disc. 

These two-a models, although not without problems 



(Smak 200C; Lasota 2001), have basic properties (recur- 
rence time, amplitude and length of some type of out- 
bursts) which compare reasonably well with observations, 
at least for some classes of CVs (for instance, the U Gem 
type dwarf novae) . The standard model as applied to CVs 
offers a background upon which more detailed scenarios 
have been constructed (with various success) to explain 
observational features such as the UV delay, supe rout- 
bursts, superhumps, standstills etc. (see Lasota P00l| for a 
review) . None of these features have (up to now) threat- 
ened the basic premise that the outbursts are triggered by 
the thermal-viscous instability in a two-a thin disc. 

In contrast, the standard DIM was never able to re- 
produce even the most basic features of any SXT. In par- 
ticular, the impossibility of obtaining anything like the 
observed exponential decays casted serious doubts on the 
validity of using the two-a assumption to represent the 
viscosity in these discs, prompting the use of more com- 
plex formulations (Cannizzo et al. 1995] ). After all, disc 
models with ab initio treatment of the viscosity are still 
some way off and the a-disc could be an oversimplifica- 
tion. However, one must ensure that the failure of models 
to reproduce observations are really due to the viscosity 
prescription, and not to the fact that some physical effects 
have not been correctly taken into account. 

Indeed, in a number of cases it has been shown that 
the disagreements between models and observations re- but very unlikely for black-holes: Lasota 200C) or the disc 



shown by Smak (1998) that the failure of models to repro- 
duce the observed UV delay largely resulted from the use 
of a fixed outer boundary condition and too small discs, 
making outside-in outbursts quite difficult to obtain. And 
of course, it was early realized that parameters such as the 
mass transfer rate from the secondary had to vary in order 
to account for the standstills of the Z Cam stars (Meyer 
& Meyer-Hofmeister 1983), not to mention the necessity 
of using a number of grid points large enough to resolve 
fronts (Cannizzo 1993). It will therefore come as no sur- 
prise that the failure of the two-a model in SXTs was more 
indicative of missing physics than inadequate assumptions 
on the viscosity. 

In the following we review the results of the previ- 
ous sections to show how the fundamental properties of 
evaporation and irradiation lead to better agreement with 
observations within the two-a model. We then critically 
examine possible failures of this model in light of detailed 
comparisons with observations. 

7.1. Fundamental properties of the model 

It has long been known that the standard disc instabil- 
ity model fails on three fundamental accounts: (a) the 
timescales of the cycles are wrong with short outbursts 
and recurrence times; (b) the decay from outburst is domi- 
nated by reflare episodes which are not observed; (c) even 
the full conversion to X-rays of the rest mass energy of 
the accreted material could not power the observed lumi- 
nosities if the quiescent disc reaches down to the com- 
pact object. All of these are unavoidable consequences 
of the model which do not depend on the detailed as- 
sumptions. In the standard DIM thermal cooling starts as 
soon as the outburst peak is reached. This proceeds on a 
short timescale since an unstable ring evolves on a thermal 
timescale. This results in (a). Problem (b) is the conse- 
quence of the high densities behind the cooling front and 
of the decreasing critical density needed to trigger an out- 
burst as the cooling front propagates inwards. Unless the 
disc can be emptied in outburst, which is difficult because 
of (a), reflares are unavoidable. Problem (c) is also due to 
the decreasing critical density: a quiescent disc must have 
densities below the critical Smax which implies extremely 
low accretion rates if the disc extends to small radii. Ei- 
ther the quiescent X-rays are not produced by accretion 
(which is possible for neutron stars: Rut ledge et al. 2000 , 
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is truncated so that the mass flow through the inner region 
(assuming it is constant) is high. This was first suggested 
in the context of dwarf novae to explain the optical to 



EUV delay (Meyer & Meyer-Hofmeister [19941 ). 

This is a particularly attractive scenario when the in- 
ner region is a hot ADAF flow (Lasota et al. 1996). An 
ADAF of varying size and accretion rate produces harder 
or softer spectra which compare well with the observed 
X-ray spectral states of SXTs. With this in mind, Menou 



et al. ( |2000| ) included disc truncation in quiescence to the 
standard DIM. Yet, their models required values of a lower 
than those usually assumed in dwarf novae and were not 
able to reproduce the short rise times observed in SXTs. In 
principle one would hope the viscosity behaves in the same 
way between CVs and LMXBsj^. Another way to get long 
recurrence times would be to have the rate of evaporation 
of the disc in quiescence almost match the mass trans- 
fer rate from the secondary (Meyer & Meyer-Hofmeister 
199£). Again this is not satisfactory: this requires fine- 



tuning of the evaporation process and/or strongly implies 
that most of the LMXB population would be cold an d sta- 
ble (Menou et al. |l999b| ; Meyer & Meyer-Hofmeister ^OOOD 
with the outbursts triggered by some unspecified mecha- 
nism (e.g. increased mass transfer from the secondary). 
This is not unlikely, but is disappointing in that most 
SXTs would then be unusual in some respect with little 
prospect for unifying explanations. 

Irradiation provides an adequate resolution of the 
three problems mentioned without having to assume a 
very low quiescent viscosity. Irradiation is an observational 
fact and simple estimates show its impact on the outburst 
cycles must be important. Irradiation is also the only ma- 
jor difference between CVs and LMXBs. In CVs the ir- 
radiation flux is at lower energies and originates from a 
(relatively) large white dwarf. In LMXBs the flux is at 
high energies and (presumably) originates from the cen- 
tral regions of the disc or from a corona. The distribution 
and amplitude of the fluxes are totally different. As such 
it was a prime candidate for the most important missing 
piece of physics in SXT models. 

In §4 we have seen how irradiation drastically changes 
the outburst. The key point is that the thermal instabil- 
ity is suppressed when Tirr ^ 10^ K. In this sense, the 
results shown in this paper depend little on the assump- 
tions made on the irradiation flux. A disc will show slow 
exponential decays for as long as a large enough irradiat- 
ing flux can keep the outer disc above hydrogen ionization 
and the exact formulation of F\„ is unimportant. The first 
major consequence is therefore to prolong the outburst by 
providing a significant additional source of heating. 

During the thermal decay, the irradiation flux deter- 
mines the location of the cooling front. In general the cool- 



ing front can only appear after most of the mass of the disc 
has already been accreted. Even if irradiation stops too 
quickly and the outburst ends prematurely, the densities 
in the disc will still be much lower than when irradiation is 
not included. This is our second major consequence since 
this makes reflares impossible and also lengthens the qui- 
escence time needed to rebuild the disc mass. 

In our models we have neglected the possible impact 
of the X-ray spectrum on irradiation heating. We have 
implicitly assumed it was soft enough to confine heating 
to a small layer at the top of the disc (Lyutyi & Sunyaev 
1976). On the other hand, hard X-rays should penetrate 
deeper and deposit energy all along the vertical layer. This 
should not significantly change the main conclusions of 
this paper for a simple reason: whatever the distribution 
of irradiation heating in the layer it will still globally heat 
the ring and displace its thermal equilibrium S'-curve to 
lower E. This can be seen in Figs. 1-2 of El-Khoury & 
Wickramasinghc ( 1999| ). Even hard X-rays will prolong 
the outburst until sufficiently low densities can be reached 
for the disc to cool. 

In an irradiated disc extending to low radii problem 
(c) remains. By including truncation in §5, we have shown 
that very long recurrence times are easily obtained. These 
result naturally from irradiation which forces the disc to 
accrete more during outburst and from evaporation which 
requires the disc to reach higher densities before an out- 
burst can be triggered again. There is no need for very low 
values of a in quiescence to reach recurrence timescales of 
tens of years. 

An adequate model of SXTs can therefore be obtained 
by including some form of irradiation in outburst and some 
form of evaporation in quiescence. Regardless of the de- 
tailed assumptions this will necessarily lead to long out- 
bursts with slow decays and long quiescence times for val- 
ues of the viscosity which are similar to those used in 
models of CVs. 

1.2. Comparison with observations 

The observed timescales for the rise (1-14 days), decay 
(8-40 days), and total outburst duration (40-300 days) 
defined by Chen et al. (1997) are well reproduced by our 



models in very general conditions i.e. with no fine tuning of 
any parameter, including a. The model also predicts recur- 
rence timescales anywhere between a year to tens of years, 
as observed. The amplitudes of the outbursts are high with 



(assuming = O.lMpeakC ), all 



^ But one should keep in mind that there is a subclass of 
dwarf novae (the WZ Sge systems) for which low values of the 
viscosity in quiescence have been invoked 



peaks at lO^-'-ss g^g g-i 
consistent with observations. Weak outbursts (i.e. with 
less irradiation) have shorter durations than strong ones 
just as noticed by Chen ct al. (1997). 

Our models, taken at face value, do not make accu- 
rate predictions for the observed outburst amplitudes in 
neutron stars and black holes because the efficiency e of 
the ADAF in quiescence is not properly calculated here. It 
would be particularly interesting to see whether quiescent 
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days 



Fig. 20. V-band optical flux (full line) compared to Afjir ~ 
eAfin (dotted line) during an outburst of the model of §5 
(Fig. ^0|). Mirr is taken as an indicator of the X-ray flux. 
Both are normalized to their peak value. The optical flux 
is halved in « 100 days while the X-ray flux is halved in 
« 60 days. 



discs with neutron star primaries have similar luminosities 
as black hole systems as briefly suggested in §3.5. In this 



case, this would support the idea that the higher quiescent 
luminosities of NS SXTs is due to thermal emission from 
the surface of the NS heated by accretion during the out- 
burst and hence proves indirectly the existence of horizons 



in BHs (Narayan et al. 1997; Menou et al. 1999c) 



A clear feature of the model which is observed in SXTs 
is the slower optical decay time. This is due to irradiation 
which, under very general conditions, can keep the cold 
outer disc hotter than usual during the propagation of the 
cooling front (Fig. Fig. |lj shows the disc behind the 
front is dominated by irradiation and almost isothermal. 
Obviously, the peak optical flux is much higher when the 
disc is irradiated than when it is not. 

There is observational evidence pointing towards a de- 
lay between the outburst rise in the optical and the rise in 
the X-rays. This has been most evident in GRO J1655-40 
where a 6 day delay was observed (Orosz et al. 1997). As 
stressed by Hameury et al. 1997 , this delay and the light- 
curve shape cannot be explained by an outside- in outburst 
in a disc extending down to the last stable orbit. The delay 
is rather the time for the thin disc to replace the hot flow in 
the inner regions. This happens on a much longer viscous 
timescale than the timescales on which fronts propagate 
through the disc. In the early rise, one would expect the 
optical to rise flrst due to the propagation of the heat front 
and the X-rays to form later when the thin disc reaches re- 
gions close to the compact object. Whether the outburst 
starts with an inside-out or an outside-in front matters 
little. 

In the present calculations, the time taken by the thin 
disc to reach its minimum radius is of the order of a few 
days (see §5.1). Although it would be desirable, we can- 



not provide better tests of the delays as this would re- 
quire a detailed model of X-ray formation and of irradi- 
ation during the early rise (which can drastically change 



the optical lightcurve). This can be done if one couples 
an ADAF model with the thin disc and consistently cal- 
culates the evolution of both, specifically the thin disc 
irradiation from the ADAF. 

To zero order the model thus satisfies all major ob- 
served properties. Inversely, what can be inferred from an 
observed lightcurve on the accretion disc ? An easily ob- 
tained measure is the peak flux of the outburst. The study 
of §6 shows this increases with mass accretion rate and disc 
size (orbital period) . The peak luminosity for neutron star 
primaries is generally smaller than for a black hole with 
a similar accretion disc. However, this might be depen- 
dent on the assumed evaporation law (here evaporation 
is weaker for a NS). Another parameter of interest, when 
available, is the recurrence time between outbursts. The 
models predicts the lowest times for systems closest to the 
stability limit i.e. with high mass transfer rates for their 
sizes (or vice- versa). Systems with long orbital periods are 
more likely to show outbursts deviating from an exponen- 
tial since the heat front may not reach the outer edge. The 
strength of irradiation is also important : systems showing 
little evidence for it would have short outbursts. The above 
assumes identical values of the viscosities. In practice, the 
effects of the different parameters are too entangled to 
reach strong conclusions from an observed lightcurve ex- 
cept, perhaps, in the most extreme cases. Consequently, 
the prospects of further constraining viscosity from the ob- 
served outburst lightcurves are dim. Dwarf novae, if only 
because of shee r num bers, are a better tool to achieve this 
goal (e.g. Smak 1999 ). X-ray observations of SXTs in qui- 
escence will also be very helpful in constraining both the 
evaporation processes and the viscosity at low tempera- 
tures, since very low quiescent viscosities should result in 
very small X-ray luminosities. 

7.3. Outstanding problems 

A quick review of observed X-ray lightcurves will reveal 
that these are much more complex than the nicely behaved 
decays we predict. Indeed, the past few years have seen 
a number of SXTs detected with lightcurves which seem 
destined to ruin a model originally based on the canoni- 
cal A0620-00 (e.g. GRO J1655-40, XTE J1550-564, XTE 
J1118-I-480...). These systems typically show consecutive 
outbursts separated by a few months with the second out- 
burst more akin to a plateau than an exponential. Both 
the plateau and short quiescence time (if the system re- 
turns to quiescence) cannot be explained by the model pre- 
sented here without additional assumptions. Possibilities 
include increased mass transfer from the secondary after 
the star was irradiated during the flrst outburst. This can 
easily lead to plateau-like lightcurves when mass transfer 
temporarily stabilizes the disc. Further complications can 
arise if the disc is large enough that it becomes warped in 
outbu rst and , for instance, sha dows the secondary (Esin 
et al. |20004 Ogilvie & Dubus |2000[ ). Winds created by 
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strong irradiation can cause significant mass loss resulting 
in quicker exponential decays (Mineshige et al. 1993| ). 

The most decently behaved SXTs also have some in- 
triguing features. For instance, A0620-00 and GRS 1124 - 
68 show glitches (in the terminology of Chen et al. 1997) 
during a decay which otherwise could very well be fitted 
by the model (see the discussion by Cannizzo 1998b). Our 
models clearly show these cannot be due to the transition 
between the viscous and thermal decays as speculated by 
King & Ritter ( |1998D . Esin et al. ( |20()0b[ ) proposed that 
the glitches could be due to increased irradiation in a 
warped disc while Cannizzo ( 200C| ) argued for evapora- 
tion in the outer disc. Both these possibilities have been 
explored and will be dealt with in a future paper. Even 
more puzzling are the multiple 'mini-outbursts' observed 
at the end of the outburst of e.g. GRO J0422-f 32. These 
bear many similarities to so me ou tbursts of WZ Sge type 
dwarf novae (e.g. Kuulkers 2000 ), suggesting a common 
mechanism. The mini-outbursts are in no way compara- 
ble to the reflares discussed above. There are no satisfying 
explanations for the mini-outbursts. 

8. Conclusion 



We have shown that, under general assumptions, a model 
of the thermal viscous instability including both the effects 
of irradiation and evaporation of the inner disc in quies- 
cence led to outburst cycles with timescales, amplitudes 
and slow decays which compare well with observations. 
This does not require any fine-tuning or modifications of 
the standard viscosity assumptions used in models of CVs. 
Even very moderate amounts of irradiation significantly 
increase the length of the outburst, its shape and the re- 
currence timescale. 

These results depend little on the detailed choices 
made in describing either the irradiation flux or the evap- 
oration law. The assumptions are most important during 
the rise and at the end of the outburst. Further compar- 
isons with observations will require these to be known 
more precisely. The next logical step will be to consis- 
tently include an ADAF model in the evaporated region 
so as to be able to calculate at each moment the X-ray 
spectrum and flux seen by the disc (and the observer). 

The wide applicability of the DIM in CVs and SXTs 
is therefore conflrmed. The importance of a careful as- 
sessment of the necessary physics cannot be more stressed 
than for the case of SXTs whose basic features can only be 
reproduced by the model when the effects of irradiation 
and evaporation are included. 

Yet, the lightcurves of SXTs show more complexity 
than the simple exponential or slow decays we predict. 
This is suggestive of still more missing physics, possibly 
variations in the mass transfer rate from the secondary 
or complex behaviour of the irradiation flux. These are 
difficult to model and most certainly vary from object to 
object. In this sense, more detailed modelling might not 



prove very rewarding. Also the prospects of using SXTs 
to constrain viscosity beyond the simple two-a model are 
poor. The flner effects of viscosity are hidden amongst the 
complex interactions between other physical effects. How- 
ever, the dynamical models show viscosity must behave in 
the same way in CVs and SXTs which, at least, supports 
a common origin. 
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Fig. 5. The outburst rise for the 
model of §4 (irradiated, with a fixed 
i?in). The upper left panel shows 
Mill and Min. (dotted line); the bot- 
tom left panel shows the V magni- 
tude. The spike at t « 1 day corre- 
sponds to the arrival at i?in of the 
inward propagating front. The rise 
is dominated by the outward front 
(inside-out outburst). Each dot cor- 
responds to one of the S and Tc pro- 
files in the right panels. The disc ex- 
pands during the outburst to trans- 
port the angular momentum of the 
material being accreted. The pro- 
files close to the peak are those of a 
steady-state disc (S oc Tc i?"'^/^). 



Fig. 6. The outburst decay for the 
model discussed in §[1.2. The 
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fully ionized disc decays viscously 
(exponential lightcurve) until t « 
70 d where Tnr{Rout) « lO"^ K. 
Thereafter the outer disc can fall on 
the cold branch of the S'-curve and 
a cooling front appears, propagat- 
ing inward only as far in as Ti„ « 
lO"* K (irradiation-controlled linear 
decay). At t « 190 d irradiation 
shuts off and the disc cools quickly. 
The dashed line is Smin for a^, with- 
out irradiation, showing the post- 
front S are very low when irradi- 
ation is included. Smin,max(a, Tin), 

and Tin- are shown for the last S and 
Tc profiles (dotted lines). 
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Fig. 7. Quiescence for the model 
discussed in ^4.6. The irradiation 



flux is negligible in quiescence and 
the overall evolution is the same as 
in the standard DIM. Mass transfer 
from the secondary is slow enough 
that matter diffuses down the disc, 
gradually increasing Mi^. The out- 
burst is triggered inside at i? w 2 x 
10^ cm when E reaches Smax (dot- 
ted line). Lower densities at the be- 
ginning of the quiescent state (due 
to the irradiation-controlled out- 
burst decay) lead to a long recur- 
rence time (about 10 years). 
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Fig. 11. The outburst rise for the 
model of §5 which includes irradi- 
ation and evaporation (see §5.1). 
Upper left panel shows Afi„ and 
Mirr (dotted line) ; bottom left panel 
shows the V magnitude. Right pan- 
els show radial profiles at the dif- 
ferent times indicated by the black 
dots on the curves. Evaporation de- 
creases during the first part of the 
rise with the thin disc extending 
to smaller radii. At i ~ 5.5 days 
the thin disc reaches the minimum 
possible inner radius of the model. 
The disc then behaves in exactly the 
same way as in Fig. ^ 
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Fig. 12. The outburst decay for the 
model discussed in The disc be- 
haves in the same way as in Fig. |^ 
until evaporation sets in at t « 

170 days (Mi„ = Movap(i?min)). 

This cuts off irradiation and the disc 
cools quickly. In contrast to Fig. |[ 
the irradiation cutoff happens be- 
fore the cooling front could propa- 
gate through most of the disc, hence 
the irradiation-controlled linear de- 
cay (i w 80 — 170 days) is less obvi- 
ous in the time profile. Ti„ (dotted 
line) is shown for the last tempera- 
ture profile. 





H (cm) 



Fig. 13. Quie scen ce for the model 
discussed in §5.2. The irradiation 
flux is negligible in quiescence. Over- 
all, the evolution is the same as in 
the standard DIM. Mass transfer 
from the secondary is slow enough 
that matter diffuses down the disc, 
gradually increasing Afjn (and low- 
ering i?in). The outburst is trig- 
gered at i? « 10^" cm when S 
reaches Smax (dotted line). Lower 
densities at the beginning of the qui- 
escent state (due to the irradiation- 
controlled outburst decay) lead to 
a long recurrence time (about 35 
years) . 



